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Abstract: 

We study the impact of electroweak radiative corrections on direct production of 
photons with high transverse momenta at hadron colliders. Analytic results for the 
weak one-loop corrections to the parton scattering reaction qq — > jg and its crossed 
variants are presented. For the high-energy region, where the corrections are strongly 
enhanced by logarithms of s/Myy, we derive simple asymptotic expressions which 
approximate the exact one-loop results with high precision. The dominant two-loop 
electroweak contributions are also calculated. Numerical results are presented for 
the LHC and the Tevatron. The corrections are negative and their size increases 
with transverse momentum. For the LHC, where transverse momenta of 2 TeV or 
more can be reached, the size of the one- and two- loop effects amounts up to —17% 
and +3%, respectively. At the Tevatron, with transverse momenta up to 400 GeV, 
the one-loop corrections do not exceed —4% and the two-loop effects are negligible. 
Finally we compare the cross sections for hadronic production of photons and Z 
bosons and find that the electroweak corrections have an important impact on their 
ratio. 
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1 Introduction 



The study of direct photon production, consisting of the QCD Compton process 
gq — > 75 and the annihilation process qq — > ~fg, has always been an important topic 
in hadron colliders physics, both theoretically PHEJ arid experimentally jH|. Since the 
photons do not fragment and can be clearly identified experimentally, direct photon 
production provides a much clearer probe of the hard-scattering dynamics than 
jet production processes. Therefore the study of large transverse momentum direct 
(prompt) photon production constitutes an important test of perturbative QCD and 
the point-like nature of quarks and gluons. Furthermore, the reaction contributes to 
background for many signals of new physics. Being embedded in the environment 
of hadronic collisions, the reaction necessarily involves hadronic physics, like parton 
distributions, and depends on the strong coupling constant. In turn, the cross 
section for direct production of photons, their tranverse momentum (px) and rapidity 
distributions can be used to gauge the parton distribution functions. Since the 
gluon distribution enters already at the leading order, the measurement of the direct 
photon production is an important means to constrain information on the gluon 
content of the proton jlj. In particular, large transverse momentum production 
provides a unique opportunity for determination of gluon densities at large x. 

Apart from the direct process, prompt photons can be also produced through 
a fragmentation process. However, most of the fragmentation contribution can be 
removed by applying an isolation criterion. The importance of the remaining contri- 
bution from fragmentation, after applying the isolation cut, is expected to decrease 
with higher px- Moreover, background processes to isolated direct photon produc- 
tion, i.e. photon production through decays of neutral mesons (71*°, rj) coming from 
jet fragmentation, are shown to be less important at large px P2|- To achieve reliable 
predictions at high px, QCD corrections in next-to-leading order |2j, are manda- 
tory. The corrections can amount to several tens of percent [3 El El, depending 
on the observable under consideration, value of px, an d details of the calculation 
such as jet definition or renormalization and factorization scales. The evaluation of 
next-to-next-to-leading order corrections involves two-loop virtual plus a variety of 
combined virtual plus real corrections and is a topic presently pursued by various 
groups [Zj. 

For the experiments at the Large Hadron Collider (LHC) a new aspect comes into 
play. The high center-of-mass energy in combination with the enormous luminosity 
will allow to explore parton-parton scattering up to energies of several TeV and 
correspondingly production of gauge bosons with transverse momenta up to 2 TeV 
or even beyond. In this region electroweak corrections from virtual weak boson ex- 
change increase strongly, with the dominant terms in L-loop approximation being 
leading logarithms of the form a L log 2L (s/M^), next-to-leading logarithms of the 
form a L log 2i_1 (s/M^), and so on. These corrections, also known as electroweak 
Sudakov logarithms, may well amount to several tens of percent. They have been 
studied in great detail for processes involving fermions in Refs. (HI E] ■ Investigations 
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on the dominant and the next-to-leading logarithmic terms are also available for 
reactions involving gauge and Higgs bosons [HI HU E2| • A recent survey of the liter- 
ature on logarithmic electroweak corrections can be found in Ref. jT3|. The impact 
of these corrections at hadron colliders has been studied in Ref. [Hj. Specifically, 
hadronic Z-boson production at large pr has been investigated in next-to-leading 
logarithmic approximation, including the two-loop terms [Tj3 22] • Numerical results 
for the complete one-loop terms have been presented in Ref. both for Z and 
photon production. 

It is the aim of this work to obtain an independent evaluation of the complete 
one-loop weak corrections to the photon production at large transverse momentum, 
and to present the full result in analytic form. At the partonic level the reactions 
qq — > 7<7, gq — > 7(7 and gq — > 75 with q = u,d,c,s or b have to be considered which 
are, however, trivially related by crossing and appropriate exchange of coupling con- 
stants. Much of the calculation proceeds in close similarity with the evaluation of 
Z production presented in Ref. [T£|. In fact, a significant part of the result can 
be obtained directly from the expressions published in Ref. [TS]. We split the cor- 
rections into an "Abelian" and a "non-Abelian" component. We present analytic 
results for the exact one-loop corrections that permit to predict separately the var- 
ious quark-helicity contributions. We also derive compact analytic expressions for 
the high-energy behaviour of the corrections. Here we include quadratic and linear 
logarithms as well as those terms that are not logarithmically enhanced at high ener- 
gies but neglect all contributions of 0(M^/s). The accuracy of this approximation 
is then discussed. 

After convolution with parton distribution functions, radiatively corrected pre- 
dictions for transverse momentum distributions of photons at hadron colliders are 
obtained. Concerning perturbative QCD, these predictions are based on the lowest 
order and thus proportional to as- To obtain realistic cross sections, higher-order 
QCD corrections would have to be included, in next-to- leading or even next-to- next- 
to-leading order. 

The paper is organized as follows: In Sect. [21 the Born approximation, our con- 
ventions and the kinematics are introduced. Sect. 01 contains a detailed description 
of the radiative corrections. Sect. 13. ll is concerned with the renormalization proce- 
dure. In Sect. 13.21 our one-loop results are given in terms of expressions that had 
been introduced in Ref. for the case of Z-boson production. Special attention is 
paid to results in the high-energy limit (Sect. 13. 3[) . which can be cast into a compact 
form. Sect. [U contains results for the dominant two-loop terms, i.e. the leading and 
next-to-leading Sudakov logarithms. In Sect. El a detailed discussion of numerical 
results is presented, both for pp and pp collisions at 14 TeV and 2 TeV, respec- 
tively. The leading and next-to-leading two-loop logarithmic terms are included in 
this numerical analysis. Sect. [Ul contains a brief summary. 
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2 Preliminaries 



The transverse momentum (pp) distribution of photons in the reaction hih 2 —>■ 7+jet 
is given by 



da hlh2 



V / dxi / dx2 6(x 1 X2-f min )f hui (x 1 ,fi 2 )f h2 j(x 2 ,n 2 )- — , (1) 

where f min = Ap^/s and y/s is the collider energy. The indices i,j denote ini- 
tial state partons (q,q,g) and fh lti (x, /i 2 ), fh 2 ,j(x,fi 2 ) are the corresponding parton 
distribution functions. The partonic cross section for the subprocess ij — > jk is 
denoted by cr* J and the sum in $1} runs over all i, j combinations corresponding to 
the subprocesses 

qq -» 7#, 9? TP: 99 -» 79, qg -» 79, 99 -» 79, 99 -» 79 (2) 

with q = u,d,c, s or 6. The Mandelstam variables for the subprocess ij — > 7/c are 
defined in the standard way 

s = (pi + pj) 2 , i = (pi - p 7 ) 2 , u = (pj - p 7 ) 2 . (3) 

The momenta p, pj, pu and p 7 of the partons are massless. In terms of Xi,x 2 ,Pt 
and the collider energy y/s we have 

s s 
s = Xix 2 s, t — — -(1 — cosd), u = — -(1 + cos6), (4) 



with cos 6 = yl — Ap^/s corresponding to the cosine of the angle between the 
momenta pi and p 7 in the partonic center-of-mass frame. 

The angular and the pr distribution for the unpolarized partonic subprocess 
ij — > 7 A; read 

d& " 1 El^l 2 (5) 



dcos# 32irNijS 



and 



da*- 2 pt 



dpT 8irNijs\t — u 



^2\M ij \ 2 + (t 



(6) 



where X) = | S po i S co i involves the sum over polarization and color as well as the 
average factor 1/4 for initial-state polarization. The factor l/N^ in ©-©5 with 
N qq = N q q = N 2 , N gq = N qg = N qg = N gq = N C (N 2 - 1), and iV c = 3, accounts for 
the initial-state colour average. 

In the following we present analytic results for the unpolarized squared matrix 
elements X)|-M^ 9 | 2 for the qq — > 7^ process. Corresponding results for the other 
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Figure 1: Tree-level Feynman diagrams for the process qq — > 7(7. 



processes in (J2j) are easily obtained by means of CP symmetry and crossing trans- 
formations [TT)j . 

The lowest-order contribution results from the tree diagrams depicted in Fig. Q 
and reads 



where Q u = 2/3, Qd = —1/3, a = e 2 /(47r) and «s = oili.^ 71 ) are the quark electric 
charges, and the electromagnetic and the strong coupling constants, respectively. 

3 One-loop corrections 

In this section, we present analytic results for the one-loop weak corrections to the 
process qq — > jg. They are closely related to those presented in Ref. (THj for the 
process qq — > Zg. In fact, the one-loop diagrams for qq — > jg, see Fig. EJ are 
analogous to the diagrams for qq — > Zg. Hence a large subset of corrections to 
<?<7 — *■ 19 can De obtained directly from the results for qq — > Zg by setting M z to 
zero. However, the renormalization counterterms need to be calculated separately 
for the qq — > jg process. 

Similarly as in the Z-boson production case [16 , we do not include electromag- 
netic corrections in the full one-loop result. The renormalized weak corrections 
are infrared finite. To obtain the results of Ref. ^j, we assumed all quarks to be 
massless, neglected diagrams involving couplings of quarks to Higgs bosons or would- 
be-Goldstone bosons and omitted quark-mixing effects. The calculation of Ref. ^Ej 
was performed at the level of matrix elements and allows for a full control over po- 
larization effects. The one-loop amplitude was split into two parts according to the 
structure of the gauge group generators in front of each term: Abelian (characteris- 
tic for Abelian theories) and non-Abelian (originating from the non-commutativity 
of weak interactions). Tensor loop integrals, appearing in the expressions for one- 
loop corrections were reduced to scalar integrals by means of the Passarino-Veltman 
technique. Assuming zero quark masses led to mass singularities of collinear nature, 
present at the level of individual integrals. However, these cancel between various 
scalar integrals, making it possible to write the final result in terms of finite combi- 
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(7) 
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Figure 2: One-loop Feynman diagrams for the process qq — > 7*7. The diagrams 
v5, v6 and b3 involve only charged weak bosons, W , whereas the other diagrams 
receive contributions from neutral and charged weak bosons, V — Z, W ± . 



nations of the integrals. Further details concerning the calculation can be found in 
Ref. PS]. 

In the following, only these contributions are presented in detail that cannot be 
directly derived from Ref. [THj by setting Mz to zero. We start discussing renor- 
malization in the on-shell scheme, in Sect. 13.11 The complete one-loop corrections 
and their asymptotic behaviour are then presented in Sect. 13.21 and Sect. 13.31 re- 
spectively. 

3.1 Renormalization 

The renormalization of the qq —>■ jg process is provided by the ^qq counterterm. 1 
This can be written, similarly to the Zqq counterterm as 




1 Samc as for the case of the qq — > Zg process the counterterm contributions associated 
with the gqq vertex and the quark propagator cancel. 
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with uj r = (1 + 7 5 )/2 and uj l = (1 - 7 5 )/2. As in Ref. HE], 



5Z, 



v=z,w± 



7A 



M 2 



(9) 



are the wave-function renormalization constants of massless chiral quarks. The 
coupling factors in @ read 



v y <?a V <• 



■ '/A 



c w 2 



v=w ± 



jVjV 
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XL 



<1\ 



2s 2 



(10) 



with the shorthands c w = cos^w and s w = sin^w for the weak mixing angle #w 
Y qx /2 = Q q — Tg A and T Q 3 A is the weak isospin of chiral quarks, i.e. T^ A = <5al/2 



and T| A = — <5al/2. The remaining counterterm contributions in (jHJ are split in an 
Abelian part, SC^ x , proportional to the tree-level coupling Q q and a non-Abelian 
part, SC^ X , proportional to the weak isospin In the following we adopt the 



on-shell renormalization scheme jTH], with = 1 — = M^jM\ and the electro- 
magnetic coupling constant defined in the Thompson limit, i.e. at the scale zero. In 
this scheme, the Abelian counterterm 5C q vanishes owing to Ward identities [TKj . 



5C. 



7A 



\z AA + S ^5Z ZA + 6 4)=0- 



The expressions for the counterterms 5Z AA , 5Z ZA and 6e 2 /e 2 can be found in 
Ref. pH]- As well known, the on-shell counterterm Se 2 /e 2 contains large logarithms 
of light-fermion masses, which are responsible for the running of a from the scale 
zero to the characteristic scale of the process. However, owing to the Ward identity, 
for the case of on-shell photon production these logarithms are cancelled by corre- 
sponding terms present in 5Z AA , see (JTTJ). This justifies our choice of a at the scale 
zero as input parameter. 

The non-Abelian counterterm reads 



1 



2 S \y C\\r 



-5Z 



ZAi 



where 



with 



5Z 



ZA 



S£(0) 



a 4c v 

47T S w 



A 



uv 



log 



1/e - 7e + log(47r) + log 



Ml 



(12) 



(13) 



(14) 



in D = 4 — 2e dimensions. Here /xr denotes the scale of dimensional regularization. 
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3.2 Result 



The complete 0(a 2 as) result for the qq — > process can be written as 



qq\2 



87T 2 aa s (iV c 2 - 1) 




Z7r y=z,w± 



(15) 



where ifn 



(t 2 + m 2 ) / {iii) is the Born contribution and 5C^ x 

r A/N, 



is the counterterm 

specified in Sect. 13. ll The remaining functions Hi (My) represent the Abelian 
(A) and non- Abelian (N) contributions resulting from the loop diagrams (see Fig. EJ) 
involving virtual weak bosons, V = Z, W , and the fermionic wave- function renor- 
malization constants Q. The above result is analogous to the one presented in 
Ref. [IE] for the qq — ► Zg process. In particular, the functions H^ N (M v ) for the 
qq — > 7<? process can be obtained from the corresponding functions for the qq — ► Zg 
process [16J by setting the mass of the external Z boson equal to zero. 2 This simple 
relation between the loop corrections to the processes qq — *> 75 and qq — > is due 
to the fact that, for Z production, the longitudinal contributions to the Z-boson 
polarization sum vanish. This is a consequence of the Ward identities presented in 
Sect. 3.6 of Ref. 



3.3 High-energy limit 

Let us now consider the high-energy region, s/Myy ^> 1, where the weak correc- 
tions are strongly enhanced by logarithms of the form log(s/Mjy). As in Ref. ^Ej 
we derive compact analytic expressions that describe the asymptotic high-energy 
behaviour of the one-loop corrections, H^ N (M v ), to next-to-next-to- leading loga- 
rithmic (NNLL) accuracy. This approximation includes double and single logarithms 
as well as those contributions that are not logarithmically enhanced at high energies. 
The NNLL expressions for Hf(M^), H^(M^) and H^(M 2 Z ), obtained by means of 
general results for the high-energy limit of one-loop integrals ^H], have the form 



Hi (My) NN = LL Re 



% K M v)^-, — + 9i ( M v0^ — + 92 {My) 



(16) 



In principle, the functions g i (My) for the qq — > jg process can be obtained from 
the corresponding functions for the qq — *> Zg process ^E] by setting the mass of the 



2 The functions H^ N (My) presented in Ref. {H^ consist of linear combinations of loop integrals 
Jo, ... , J14. We observe that the integrals J2, J& and J10 are logarithmically singular for Mz — ► 0. 
However the corresponding coefficients are proportional to M|. Thus the contributions of such 
loop integrals vanish for Mz ~ * 0. 
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external Z boson equal to zero. However, this must be done with care since the 
Z-mass dependence of the gf^{M v ) functions in Ref. [THj is not shown explicitly. 
For the case of photon production we find that the functions g^(M v ) with i = 0,1,2 
and gf(Mw) with j = 1,2 are the same as in Ref. JH]. Instead, for the non-Abelian 



function g${M$y) 



we obtain 



g^{Ml) = 2 A uv 



log 



Ml 




+ log 2 



t 



log 2 



M 2 W . 



-log 2 



— u 
Mr. 



l og 2 ( J j + log 2 | 



2tt 2 



(17) 



Here the logarithmic terms are the same as in the Z-production case whereas the 
constants are different. In particular, the coefficient of the 7r 2 term is different and 
no term of the form is present in (fTTjl . Logarithms with negative arguments 

in (fT7|) are defined through the usual is prescription, f — > f + ie for f = s,t, u. 
Assuming, as in Ref. ^H], that one of the Mandelstam invariants is positive and the 
other two are negative, for the real part of ()17|) we obtain 



Reg^M, 



w) 




log 2 



\u\ 



M 2 W , 



This result applies directly to the qq — > jg process, where s > 0,t < and u < 0. 
In this case 0(—s) = 0. The 9 term becomes non-vanishing when ()18|) is translated, 
by means of permutations of the Mandelstam invariants (crossing transformations), 
to the processes in (J2J) that involve gluons in the initial state. 



4 Two-loop corrections 



In the TeV energy region, the logarithmic electroweak corrections are very large. 
The one-loop contributions can amount to tens of percent and also the two-loop 
logarithmic terms can reach the level of several percent, see e.g. Ref. ^E]. Such 
higher-order contributions are thus mandatory for precise theoretical predictions. 
In the following, we present the two-loop corrections to the process qq — > jg to 
next-to-leading logarithmic (NLL) accuracy. For a discussion of the calculation we 
refer to Ref. [T^], where the same class of corrections has been computed for the 
process qq — > Zg. 

The unpolarized squared matrix element for qq — > ^g, including NLL terms up 
to the two-loop level, has the general form 



37r 2 aa s (iV c 2 



1 



t 2 + u 2 
tu 



(19) 
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The Born contribution reads 



A® = 2Ql 



(20) 



At one loop, the NLL part consists of double- and single- logarithmic terms and reads 



NLL 



E Q q 

A=L,R 



(CT - Qt) ( L * - 3L t) + ( L ? + L l ~ ^ 



(21) 



WW 



Here we used the shorthand L~ = log fc (|f|/M^) for the logarithms and C| 
Yq : [/(^ c w) + ^q\/ s w are the eigenvalues of the electroweak Casimir operator for 
quarks, with C qh = 3/4 and C QR = 0. Only weak corrections are included in (pHf . 
This result is obtained in the Mz = My/ approximation and is consistent with 
the leading- and next-to-leading logarithmic part of the one-loop asymptotic ex- 
pressions 3 presented in Sect. 01 At two loops, for the complete electroweak NLL 
corrections 4 we obtain 



NLL 



A=L,R I * \ s w , 



6L, 




f5^(Lf + L 4 -Lf 



<b 2 



(23) 



where b x = — 41/(6c^) and 6 2 = 19/(6s^) are the one- loop /3-function coefficients 
associated with the U(l) and SU(2) couplings, respectively. Eq. (|2*3Jl has been 
derived from the general results of Refs. ^2j 120] for leading- and next-to-leading 
electroweak two-loop logarithms. 



5 Numerical results 

In this section, we present numerical predictions for photon production at high 
transverse momenta, both at the partonic and hadronic level. We also compare the 
cross sections for hadronic photon and Z-boson production and study the impact of 
the corrections on their ratio. The lowest order (LO) and the next-to-leading-order 
(NLO) predictions for photon production result from (JJJ) and (|T3|l. respectively. 

3 The correspondence between (I21|) and the one-loop expressions of Sect. is easily seen by 
means of the relation 

c c q :-Ql= E ■ ( 22 ) 

4 At the two-loop level the purely weak corrections cannot be isolated from the complete elec- 
troweak corrections in a gauge-invariant way. One must thus consider the combination of weak 
and electromagnetic virtual corrections. The latter are regularized by means of a fictitious photon 
mass A = Mw- This approach is discussed in detail in Ref. 
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The next-to-leading-logarithmic 5 (NLL) approximation, at one loop, corresponds to 
the contributions (|2()jl and 1)2 lj) . The next-to- next-to-leading-logarithmic (NNLL) 
predictions, also at one loop, are based on the asymptotic expressions (ITSJl - lfTSJl 
combined with the exact result ff 13f) for the counterterm. Our best predictions, to 
next-to-next-to- leading-order (NNLO), include the exact NLO contributions com- 
bined with the leading and next-to-leading two- loop terms ()23|) . 

The hadronic cross sections are obtained using LO MRST2001 parton distribu- 
tion functions (PDFs) [25- We choose p\ as the factorization scale and, similarly 
as the scale at which the running strong coupling constant is evaluated. We also 
adopt the value ats(M§) = 0.13 and use the one- loop running expression for as(yti 2 ), 
in accordance with the LO PDF extraction method of the MRST collaboration. 
We use the following values of parameters [22]: « = 1/137.0, Mz = 91.19 GeV, 
M w = 80.39 GeV and s 2 w = 1 - c 2 w = 1 - M^/Mf . 

We begin by investigating NLO, NLL and NNLL relative corrections to the 
partonic (unpolarized) differential cross section o\a l] /d cos 9 [see ©]. To this end 
we define 

&j da^ LO /dcosfl _ 

'^NLO/LO ~~ j-ii /j n v > 

dCTLo/dCOSD 1 

and similarly T^nll/nlo anc ^ ^nnll/nlo- These ratios, calculated at cos 8 = 0, 
are displayed as a function of \f§ in Fig. El We consider four processes: uu^^g 
(Fig.Et), dd^g (Fig. Eh), gu^u (Fig.Eb), gd^d (Fig. Eli). The NLO correc- 
tions are negative, their absolute size growing with the partonic energy At 
^/I = 4TeV the corrections reach —15% for the uu and gu channels, —23% for the 
dd, and —28% for the gd channel. The NLO correction is well approximated by 
the NLL result: for \/I > 200 GeV the approximation is accurate up to 2% in all 
channels. Even better an approximation of the one-loop weak correction is provided 
by the NNLL result. In uu and dd channels the NNLL approximation differs from 
the NLO result by less than 1% in the entire region under consideration, whereas for 
the gu and gd channels the quality of the approximation is at the level of (or better 
than) one permille. The absolute LO and NLO cross section and the TZ^ ratios for 
specific values of cos#, \/I and different collision channels are listed in Tabled 

The transverse momentum distribution dcx/dpT at the LHC is shown in Fig. 0] 
We display separately the absolute values of the LO, NLO, NLL and NNLL differ- 
ential cross sections (Fig. H^,) and the relative correction wrt. the LO result for the 
NLO, NLL and NNLL distributions (Fig. HJd). The relative correction wrt. the LO 
is now defined as 

^had _ dq N Lo/dpT 1 , _v 

/^NLO/LO- d(jLo/dpT 1 ^) 

for the NLO case, and similarly for the NLL and the NNLL cross sections. The qual- 
ity of the NLL and NNLL high-energy approximations of the NLO result is shown in 
more detail in Fig. Ht. The contribution provided by the NLO correction is negative 



At the NLL level, angular-dependent logarithms are treated as in Ref. |15| . 
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Figure 3: Relative one-loop corrections to the partonic differential cross sections 
da lJ /d cos 9 at cos# = for (a) uu channel, (b) dd channel, (c) gu channel, (d) gd 
channel. The solid, dashed and dot-dashed lines denote the modulus of the 71 ratios, 
as defined in the text, for the full NLO cross section, the NLL approximation and 
the NNLL approximation of the one-loop cross section, respectively. 



and increases in size with pp- It ranges from —6% at pr = 500 GeV up to —17% at 
Pt = 2 TeV. From Fig. |3Jd and Fig. we conclude that the NLL approximation 
works well, differing from the full NLO prediction by about 3 permille at low p^ arid 
by less than 1 permille at p-r ~ 2 TeV. The quality of the NNLL approximation is 
very good, at the level of accuracy of 10~ 3 or better in the entire px range. In Fig. |3] 
we show the relative size of the corrections in the NLO approximation (^nlq/lo' 
solid line), and in the approximation which includes the next-to-leading logarithmic 
two- loop terms (T^nnxo/lo' dotted line). These additional two-loop terms are posi- 
tive, their size increasing with p^. At pt = 2 TeV they amount to 3%, yielding the 
total (i.e. together with the NLO) correction to the LO of —14%. 

To underline the relevance of these effects, in Fig. |B] we present the relative 
NLO and NNLO corrections for the cross section, integrated over p^ starting from 
Pt = Px ut ! as a function of p T ut . This is compared with the statistical error, defined 
as Ac^at/cr = 1/vN with N = £ x a LQ . We assume a total integrated luminosity 
C = 300 ftT 1 for the LHC [23J. It is clear from Fig. that the size of the one-loop 
(two-loop) corrections is much bigger than (comparable to) the statistical error. 
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3.7343 




1000 


0.5 


6.1285 x lO" 1 


5.8616 x 10- 


1 


-43.560 


8.8071 


1.9408 




2000 


0.5 


1.4080 x lO" 1 


1.2881 x 10" 


1 


-85.188 


7.6569 


0.8512 


dd 


500 


0.0 


3.9533 x 10- 1 


3.8291 x 10" 


1 


-31.427 


14.646 


6.7828 




1000 


0.0 


9.0276 x 10~ 2 


8.3494 x 10" 


2 


-75.121 


11.381 


3.6201 




2000 


0.0 


2.0771 x 10- 2 


1.7810 x 10- 


2 


-142.55 


9.5167 


1.6553 


dd 


500 


0.5 


6.7211 x lO" 1 


6.5344 x 10" 


1 


-27.769 


19.286 


5.6758 




1000 


0.5 


1.5321 x lO" 1 


1.4266 x 10" 


1 


-68.882 


16.786 


3.0445 




2000 


0.5 


3.5200 x 10~ 2 


3.0501 x 10" 


2 


-133.48 


15.642 


1.3809 


gu 


500 


0.0 


7.4125 x lO" 1 


7.2569 x 10- 


1 


-20.985 


-0.3584 


-0.5405 




1000 


0.0 


1.6927 x lO" 1 


1.5994 x 10" 


1 


-55.125 


-0.5161 


-0.3435 




2000 


0.0 


3.8945 x 10~ 2 


3.4909 x 10" 


2 


-103.63 


-0.7276 


-0.1649 


gu 


500 


0.5 


6.3010 x lO" 1 


6.1387 x 10" 


1 


-25.756 


-2.2604 


-0.2459 




1000 


0.5 


1.4364 x lO" 1 


1.3469 x 10" 


1 


-62.295 


-2.6157 


-0.1587 




2000 


0.5 


3.3000 x 10~ 2 


2.9263 x 10- 


2 


-113.22 


-3.0577 


-0.0747 


gd 


500 


0.0 


1.8531 x lO" 1 


1.7880 x 10" 


1 


-35.167 


-12.356 


-0.3746 




1000 


0.0 


4.2317 x 10" 2 


3.8365 x 10" 


2 


-93.383 


-13.653 


-0.4196 




2000 


0.0 


9.7363 x 10~ 3 


8.0397 x 10" 


3 


-174.24 


-15.298 


-0.2382 


gd 


500 


0.5 


1.5753 x 10- 1 


1.4944 x 10- 


1 


-51.298 


-18.358 


-0.0453 




1000 


0.5 


3.5909 x 10~ 2 


3.1686 x 10" 


2 


-117.62 


-20.363 


-0.1781 




2000 


0.5 


8.2500 x 10~ 3 


6.5458 x 10" 


3 


-206.57 


-23.109 


-0.1086 



Table 1: Absolute value of the LO and NLO partonic differential cross section 
d<7 l Vdcos(9 and the ratios 7^ LO/LO , "^nll/nlo and "^nnll/nlo in permille, for the 
partonic processes uu — > jg, dd — ■> jg, gu — * 7« and gd — > jd. The running strong 
coupling a s (fi 2 ) is taken at the scale /i 2 = p\ = (1 — cos 2 9)s/ A. 
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Figure 4: Transverse momentum distribution for pp^^j at \/s = 14TeV. (a) 
LO (solid), NLO (dashed), NLL (dotted) and NNLL (dot-dashed) predictions, (b) 
Relative NLO (solid), NLL (dotted) and NNLL (dot-dashed) weak correction wrt. 
the LO distribution, (c) NLL (dotted) and NNLL (dot-dashed) approximations 
compared to the full NLO result. 
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Figure 5: Relative NLO (solid) and NNLO (dotted) corrections to the px distribu- 
tion for VV^13 at \fs — 14 TeV. 




_q 20 ' I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 

200 400 600 800 1000 1200 1400 1600 1800 2000 

vT [GeV] 

Figure 6: Relative NLO (solid) and NNLO (dotted) corrections wrt. the LO predic- 
tion and statistical error (shaded area) for the unpolarized integrated cross section 
for VV~^13 & t = 14 TeV as a function of p™*. 
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Figure 7: Ratio of the transverse momentum distributions for the processes pp— >7j 
and pp^Zj at y/s = 14TeV: LO (solid), NLO(dashed) and NNLO (dot-dashed) 
predictions. 

In Fig. [3 we plot the ratio of the px distributions of photons and Z bosons. The 
latter is computed within the MS scheme using the same input parameters as in 
Ref. PE]- Such ratio is expected to be less sensitive to theoretical errors than the 
distributions themselves, since many uncertainties such as the scale at which as is 
calculated or the choice of PDFs cancel to a large extend in the ratio. Moreover, due 
to a similar cancellation mechanism, the ratio should remain stable against QCD 
corrections. From Fig. [7|we observe that the weak corrections modify the production 
ratio considerably. The effect is the strongest at high px- In this region, the LO 
photon cross section is smaller than the cross section for Z boson production by 
about 25%. The relatively large NLO corrections for Z production, as compared to 
7 production, cause the full NLO production rates to become equal at the highest 
Pt considered here, i.e. px ~ 2 TeV. The two-loop corrections modify the ratio and 
lead to a few percent decrease at high px- 

We perform a similar analysis for direct photon production at high transverse 
momentum at the Tevatron. In Fig. |H] we show the transverse momentum distribu- 
tion (Fig. |H^i), the relative size of the corrections (Fig. IHJd) and the quality of the 
one-loop NLL and NNLL approximations (Fig. |Ht). The effects of weak corrections 
are generally much smaller for the case of the Tevatron than the LHC, with the NLO 
corrections not exceeding —4% at the highest px considered, i.e. at px ~ 400 GeV. In 
Fig. E3 the relative size of the corrections to the integrated cross section is compared 
with the statistical error expected for an integrated luminosity C = 11 fb _1 |24J. At 
the energies of the Fermilab collider, the NLO weak correction is of the order of the 
statistical error and we conclude it should be taken into account when considering 
precision measurements. The two-loop terms turn out to be negligible. The ratio of 
the pt distributions of photons and Z bosons is shown in Fig. Since the weak 
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corrections to the Z and 7 production at the Tevatron are moderate, their effect 
on the ratio is fairly small and stays within a few percent range for all values of pt 
considered here. 

Our numerical results, both for the case of the LHC and the Tevatron, are in 
agreement with the results presented in figures of Ref . [T7] . 

6 Summary 

In this work we present one-loop weak corrections to the direct production of pho- 
tons with large transverse momenta at hadron colliders. Analytical results are given 
for the parton subprocess qq — * 7$ and its crossed versions. We also present ap- 
proximate expressions valid in the high-energy region, s ^> M$y, where the weak 
corrections are enhanced by logarithms of s/Myy- The complete high-energy approx- 
imation discussed in this paper includes all large logarithms as well as those terms 
that are not logarithmically enhanced. This approximation is in very good agree- 
ment with the complete one- loop result. We also calculate the two-loop electroweak 
corrections in the next-to-leading logarithmic approximation. 

The corrections are then evaluated numerically for proton-antiproton collisions at 
2 TeV (Tevatron) and proton-proton collisions at 14 TeV (LHC) in the region of large 
tranverse momentum (pr)- The corrections are negative and their size increases with 
Pt. At the Tevatron, transverse momenta up to 400 GeV will be explored and the 
weak corrections may reach up to —4%. At the LHC, transverse momenta of 2 TeV 
or more are within the reach. In this region the one- loop corrections are significant, 
about —17%, and even the dominant two-loop logarithmic terms must be included 
in precise predictions. Finally we compare the cross sections for photon and Z-boson 
production and find that the electroweak corrections have a considerable impact on 
their ratio at the LHC. 
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Figure 8: Transverse momentum distribution for pp-^'yj at ^/s = 2TeV. (a) LO 
(solid), NLO (dashed), NLL (dotted) and NNLL (dot-dashed) predictions, (b) Rela- 
tive NLO (solid), NLL (dotted) and NNLL (dot-dashed) weak correction wrt. the LO 
distribution, (c) NLL (dotted) and NNLL (dot-dashed) approximations compared 
to the full NLO result. 
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Figure 9: Relative NLO (solid) and NNLO (dotted) corrections wrt. the LO and 
statistical error (shaded area) for the unpolarized integrated cross section for pp— >^yj 
at v/i = 2 TeV as a function of p™*. 
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Figure 10: Ratio of the transverse momentum distributions for the processes pp—^^j 
and pp^Zj at = 2 TeV: LO (solid), NLO(dashed) and NNLO (dot-dashed) 
predictions. 
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